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We study gravitational lensing by compact objects in gravity theories that can be written in a 
Post-Post-Newtonian (PPN) framework: i.e., the metric is static and spherically symmetric, and can 
be written as a Taylor series in m./r, where m. is the gravitational radius of the compact object. 
Working invariantly, we compute corrections to standard weak-deflection lensing observables at first 
and second order in the perturbation parameter s = "&,/"&e, where 1?. is the angular gravitational 
radius and -Qe is the angular Einstein ring radius of the lens. We show that the first-order corrections 
to the total magnification and centroid position vanish universally for gravity theories that can be 
written in the PPN framework. This arises from some surprising, fundamental relations among 
the lensing observables in PPN gravity models. We derive these relations for the image positions, 
magnifications, and time delays. A deep consequence is that any violation of the universal relations 
would signal the need for a gravity model outside the PPN framework (provided that some basic 
assumptions hold). In practical terms, the relations will guide observational programs to test general 
relativity, modified gravity theories, and possibly the Cosmic Censorship conjecture. We use the new 
relations to identify lensing observables that are accessible to current or near-future technology, and 
to find combinations of observables that are most useful for probing the spacetime metric. We give 
explicit applications to the Galactic black hole, microlensing, and the binary pulsar J0737— 3039. 
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I. INTRODUCTION 



Gravitational lensing is now a central field of astronomy with wide-ranging applications that relate to extra-solar 
planets, dark matter substructures, and cosmological parameters (including dark energy) 00, 0j- Theoretical studies 
have also explored lensing by black holes within the context of general relativity 0-0], braneworld gravity p^|-[T^|. 
and string theory [20L |21| . Most black hole lensing studies have focused on "relativistic" images corresponding to 
light rays that loop around a black hole, which probe gravity in the strong-deflection limit, but are extremely difficult 
to detect observationally 0, 0] . We are proposing and evaluating new possibilities for using gravitational lensing by 
compact objects (including black holes) to test various theories of gravity using current or near-future technology. 

In Paper I of this series |2l| , we introduced an analytic framework for studying gravitational lensing by a compact 
deflector with mass M., in which the lensing scenario satisfies three basic assumptions: 

[Al] The gravitational lens is compact, static, and spherically symmetric, with an asymptotically flat spacetime 
geometry sufficiently far from the lens [3j| • The spacetime is vacuum outside the lens and flat in the absence of 
the lens. 

[A2] The observer and source lie in the asymptotically flat regime of the spacetime. 

[A3] The light ray's distance of closest approach r$ and impact parameter b both lie well outside the gravitational 
radius m. = GM,/c 2 , namely, m,/ro <C 1 and m./6 <C 1. The bending angle can then be expressed as a series 
expansion in m./6, as follows: 




The coefficients Ai are independent of m./fe, but may include other fixed parameters of the spacetime. Since 
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b and m. are invariants of the light ray, JIJ is independent of coordinates. Note that the subscript of Aj 
conveniently indicates that the component is affiliated with a term of order i in m m /b. 

We applied the lcnsing framework to gravity theories that can be written in a Post-Post-Newtonian (PPN) expansion, 
meaning that the metric can be written as a Taylor series in m./r. Working invariantly, we computed the observable 
properties of the primary (positive-parity) and secondary (negative-parity) lensed images. Assuming that [A3] holds, 
we wrote the lensing observables as series expansions in e = $»/$,e, where is the angular gravitational radius and 
is the angular Einstein ring radius. For gravity theories that agree with general relativity in the weak-deflection 
limit (Ai = 4 in eqn.^l, the zeroth-order terms in the e expansions give the familiar results for lensing by a point mass 
in the weak-deflection limit of general relativity. The higher order terms give corrections to the lensing observables, 
which differ for different gravity theories. We studied general third-order PPN models, which allowed us to compute 
the weak-deflection lensing results plus the first two correction terms (order e and e 2 ). 

In Paper I we found the interesting result that the first-order corrections to two lensing observables — the total 
magnification and the magnification-weighted centroid position — vanish in PPN gravity models that agree with 
general relativity in the weak-deflection limit (i.e., A\ = 4). More generally, we found that these corrections depend 
on A\ — 4, suggesting that they nearly vanish in gravity theories that agree only approximately with general relativity 
in the weak-deflection regime (i.e., Ai w 4). Note that existing observations constrain this parameter to Ai — 
3.99966 ±0.00090 IH. 

Working more generally than in Paper I, we have now discovered the surprising result that the first-order corrections 
to the total magnification and centroid in fact vanish exactly in all PPN models. This depends on precise cancellations 
that are striking because the PPN framework covers quite a broad range of gravity theories. Understanding why the 
cancellations occur has led us to identify some new fundamental relations between lensing observables in the PPN 
framework. Some of the relations are universal for PPN models, in the sense that they hold for all values of the 
invariant PPN parameters of the light bending angle. 

These new lensing relations will play key roles in planning observing missions to test theories of gravity. First, 
the relations allow us to determine which observables are most accessible to current or near-future instrumentation. 
We give explicit applications to the Galactic black hole, Galactic microlensing, and the binary pulsar J0737-3039. 
Second, the lensing relations help us identify combinations of lensing observables that are most useful for probing the 
spacetime metric by constraining invariant PPN parameters. One of the measurable parameters is connected to the 
existence of naked singularities in certain gravity models (see Paper I), so we may have an observational test of the 
Cosmic Censorship conjecture [2i|. Third, the universal relations provide a powerful means to test the entire PPN 
framework. Any violation of relations that are universal in the PPN framework would suggest that a fundamentally 
different theory of gravity is needed (provided that assumptions [Al]-[A3] hold). 

In this paper we present the new lensing relations and use them to assess prospects for using lensing observations 
to test PPN gravity theories. Section HT1 reviews the third-order PPN lensing framework. Section ITTT1 derives the new 
relations between the image positions, magnifications, and time delays, and discusses their conceptual implications. 
Section llVI considers applications of the relations to various astrophysical settings. 



II. LENSING IN THE PPN FRAMEWORK 



In this section we review our results for lensing in the PPN framework. See Paper I for the complete analysis. 
Consider a compact body of mass M, , perhaps a black hole or neutron star, that is described by a geometric theory 
of gravity. By assumptions [Al]-[A3], it suffices to analyze an equatorial metric of the form 

ds 2 = -A(r) dt 2 + B(r) dr 2 + r 2 dip 2 . (2) 

We study metrics whose coefficients can be expressed in third-order PPN expansions, 

A(r) = 1 + 2a, (|) + 2a, + 2a, (|j + ... , (3) 

*M -!-«,(* + 4 fe (*) 2 - 8t3 (*) 3 + ..., (4) 



where <f> is the three-dimensional Newtonian potential with 



m. 

r 



(5) 
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FIG. 1: Schematic diagram of the lensing geometry. Standard quantities are identified: B is the angular position of the 
unlensed source; ■& is the angular position of an image; a is the bending angle; and cLl, ds, and cLls are angular diameter 
distances between the observer, lens, and source. The impact parameter b is an invariant of the light ray and is related to the 
angular image position by ■& = sin _1 (fo/di,). 



Section III.C of Paper I derives the light bending angle for this metric, 
angle takes the form of with coefficients 

Ax = 2(oi + 6i), 

A 2 = r2a\ -a 2 + a x bi - j + 6 2 J tt , 



A, = - 



35a? + 15of6i - 3a a (l0a 2 + b\ - 4b 2 ) + 6a 3 + b\ - 6a 2 &i - 4b x b 2 



The invariant expression for the bending 

(6) 
(7) 

(8) 



The coordinate independent quantities Ai will be called the invariant PPN parameters of the light bending angle. 
For reference, the Schwarzschild metric in general relativity has a\ = b\ = b 2 = b 3 = 1 and a 2 — a% — 0, and hence 
A x = 4, A 2 = 15tt/4, and A 3 = 128/3. 

Fig. d displays the gravitational lensing scenario. Elementary trigonometry establishes the relationship (see @) 



tan ,6 = tan?9 — D (tan$ + tan(d — $)) , 



(9) 



where D — d^s/ds- This equation agrees well with the full relativistic formalism for light propagation [9j|, so we 
take it as the general form of the lens equation. The equation has two primary (weak-deflection) solutions 34]: one 
corresponding to an image on the same side of the lens as the source; and one on the opposite side. By convention, 
angles describing image positions are taken to be positive. This forces the source's angular position to have different 
signs: B is positive when we are studying an image on the same side of the lens as the source (as depicted in Fig.0; 
while B is negative when we are studying an image on the opposite of the lens from the source. 
Following Paper I, we convert to scaled variables 



B_ 

0B 



1? 



T 

TE 



1K_ 



2e. 

AD 



(10) 



where 1), — tan 1 (m,/dr,) is the angle subtended by the gravitational radius, and r is the lensing time delay. The 
natural angular scale is given by the angular Einstein ring radius, 



<jGM.d LS 

c 2 dhds 



(11) 



while the natural time scale is 



(12) 
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We then assume that solutions of the lens equation can be written as a series of the form 



= co 



h e + 6 2 e 2 + O [ef , 



(13) 



where 9q represents the image position in the weak-deflection limit, while 9± and 9 2 give the first- and second-order 
correction terms. With these substitutions, the lens equation becomes 



= D 



-4(3 + 46>o 



D 



-A 2 + (A x + 49 2 ) 9x 



D 
W 



A{ - 3A 3 + 12 A{D9q - Ax(56D 2 9% + 30( - W a 9 2 ) 



+MD 2 9l{9l - p 3 ) + 6A 2 9x + 126»^ 2 



0(ef 



(14) 



We solve for 60, di, and 9 2 by finding the values that make each term of the lens equation vanish. From the vanishing 
of the first term we obtain the relation 



(15) 



This is the generalization of the familiar weak-deflection lens equation for a point mass to Ai 7^ 4. Its solution is 

6 a = l(VAi+P 2 +(3) . (16) 

Then requiring that the second and third terms in (|14|) vanish yields the correction terms 

A2 



A x + W 2 
1 



Ax(-3A 2 2 + 3AxA 3 - Aj(D 2 - 1) 



(17) 
(18) 



3 o (Ax + 402)3 

+ 4^-6A 2 2 + 6AxA 3 + Af(D - 2)(D - l)Jflg 
+ 8(6A 3 + A?(2 + L>(llZ?-12)))6i4 

+ QAA\D{AD - 3) 0% + 128AxD 2 0% 



Notice that the O (e) 2 term in the lens equation (|14|) does not explicitly involve /3; so 9x depends on the source position 
only implicitly through 9q. By contrast, the O (e) term in the lens equation does involve (3. In writing the expression 
for 9 2 , we have found it convenient to substitute for (3 using i|15|) ■ Then the two correction terms 9x and 9 2 are both 
written only in terms of £?o- The source position dependence could be made explicit by substituting for 9q using i|16fl . 
The signed magnification \x of a lensed image at angular position is given by 



sinS(tf) dB{d) 



sin$ 



di) 



(19) 



After taking the derivative, we change to our scaled angular variables from (|1(J|I and (|13|l . and substitute for 9x and 
9 2 using lfT7|l - ((T%|l . This yields a series expansion for the magnification, 



H = Ho + \x\£ + ii 2 e 2 + 0(e) , 



(20) 



where (see eqs. 77-79 of Paper I) 



(21) 
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Mi = 



/'•2 



-A\D 2 + 8 A 2 (&A 3 + A?(2 + 6D - 9L> 2 )) 6>g 



3(Ai-4^)(^! + 4flg)8 
- 32fl8^-12AiA 3 + A^(L>(17 J D-12)-4) ,i ^ 



(22) 
(23) 



128 [6A 3 +A\{2 + QD- 9L> 2 )) 0g - 256 A 2 !) 2 ^ 



We have again substituted for (3 using (fl5|l . 

The time delay (relative to an undeflected light ray) was derived in Section V.C of Paper I. It can be written as a 
series of the form 



T = T 



n e 



o(ey 



where 



Tl 



1 

2 

M 

40 o 



ai +(3 2 



A! (d L 0*$ 



2„q2 
E 



(24) 

(25) 
(26) 



It is possible to derive the second-order correction to the time delay, but we have found that to be less important 
than the second-order corrections to the position and magnification. 

Remark. In Paper I, we gave the relation between (3 and da only for the case A\ =4. As a result, there are 
minor changes in the expressions for 62 and {12 in Paper I when A\ 7^ 4. Equations (|16fl . I|18|) . an d (|23|l give the 
correct expressions for the general case in which A\ can take on any value. Note, however, that observationally 
Ax = 3.99966 ± 0.00090 fM. 



III. NEW RELATIONS BETWEEN LENSING OBSERVABLES 



In this section we uncover several new relations between the perturbation coefficients of the fundamental lensing 
observables, namely the positions, magnifications, and time delays of the two primary images. Some of the relations 
are universal among PPN models in the sense that they hold for any values of the PPN parameters. Others hold for 
all source positions. 



1. Position relations 



Starting from l|16|l and recalling that the positive- and negative-parity images correspond to (3 > and /3 < 0, 
respectively, we can write the positions of the two images in the weak-deflection limit as 

6t = UV^+0 5 ±\l3\). (27) 



2 

We can immediately identify two interesting relations 



n - <>o = m, (28) 

0t Oo = 4 1 • ( 29 ) 



Equation l|28|l represents our first universal relation. This relation is familiar from standard weak-deflection lensing 
in general relativity (see p. 189 of 0), but now we see that it holds for all PPN models, regardless of whether they 
agree with general relativity in the weak-deflection limit. The second equation is our first example of a relation that 
is independent of the source position. 
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Next, combining l|T7 |l -l(T5 |l and (J2H yields 



n 



n 



9+ 



The first-order relation 



for 0+ 



A 1 y/A l +0P 



3AI 



6A 2 2 - 6AiA 3 - A\{2 - 3£> 2 ) 



(30) 
(31) 
(32) 



6 1 is another source-independent relation, while the one for 0^ — 1 is source 



dependent. Both first-order relations depend on the sign of A2. The second-order relation is independent of the sign 
of A2. The dependence on the sign of A2 is interesting because in Paper I we found that in certain gravity theories 
this sign is connected to the occurrence of naked singularities. 



2. Magnification relations 

We take the magnification terms l(2*T|) - (f2*3f and write 0o in terms of [3 using (|2"7jl . This yields 
± 1 M + 2/3 2 

^ = _ 4(A 1+ /32)3/2 ' (34) 

± _ , 9A 2 + 2(A 1 + P 2 ){~6A 3 + 2A\D 2 j3 2 + A\[-2 + 3£>(-2 + 3£>)]} 
/i2 ~ 241/31(^ + ^)5/2 • W 

We can then identify three universal magnification relations: 

Met + Mo = !> Mi" - Ih =°. M2 + M2 = 0. (36) 

Recall that the sign of the magnification indicates the parity of an image, so actually gives the image brightness. 
The zeroth-order relation can be rewritten as — = 1. In other words, in the PPN framework the difference 
between the fluxes of the images (at zeroth order) always equals the flux of the source in the absence of lensing. 

The first-order magnification relation arises because the first-order correction term is the same for both images, 
but the actual magnifications have opposite signs. If A2 is positive, then /ij < 0. This makes the positive-parity 
image less positive, or fainter; but it makes the negative-parity image more negative, or brighter. (If A2 is negative, 
the opposite occurs.) Consequently, the magnifications of the positive- and negative-parity images are shifted by the 
same amount but in the opposite sense. 

In addition, combining (|16|) - (|17|l with (|22ll 1231) yields another universal relation, 

( M +0+ + Mo -0n + (mM + + Mr^) = o. (37) 

This relation will be useful when we analyze the centroid position below. 



3. Total magnification and centroid 

If the two images are not separately resolved, the main observables are the total magnification and magnification- 
weighted centroid position (e.g., [2J|)- Applying the universal magnification relations H36|l to the total magnification 
yields 

Mtot = Im + I + Im~I, 

= (Mo - Mo) + (Mi"-Mr) e + (m^-M^)^ + C(e) 3 : 

= (2^-1) + + 0{ef. (38) 

There is no first-order correction to the total magnification for all gravity theories within our PPN framework. This 
result depends on a precise cancellation between /j,J and /i^, so it is striking that it is universal for PPN models. 
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An important implication is that the total magnification would have to be measured much more precisely than some 
other observables to find corrections to the weak-deflection limit (see Sec. II VI for more discussion). 
The magnification-weighted centroid position is defined by 

_ e+\»+\-e-\»-\ _ e+j^ + &jr 
ccnt " + ~ m+-m- ■ 1 J 

Writing ± and fx in terms of their series expansions, and using the magnification relations (|36|l . yields 

e cc „t = e + eie + e 2 c 2 + o{e)\ (40) 



where 



6o = +fl °">° , (41) 

Mo - Mo 

Qi = Mjgr + MiA + ±Jh6l +mX ; (42) 

Mo ~ Mo" 



q 2 _ V^O T- A/^0 ^2 _ ^0 _ ^0 AMI V"l T w l ) J ^0 "2 ~ H'Q "2 J 



9 + + #o )(M f M 2 - Mo M2 ) + (Mo -Mo)^ ^ + h ) + Mo^a" + Mo 

(Mo" "Mo ) 2 

Writing (|41|> and l|43|l in terms of {3 yields the zeroth-order centroid position and the second-order correction to be 



|/3|{9A 2 - 2(Aq + /? 2 )[6A 3 + 8A^ 2 /3 4 - 6A 2 J/? 2 (2 - 3£>) + A\{2 - D 2 )]} 
2 " 6(A 1 +/3 2 )(A 1 +2/3 2 ) 2 ■ 1 j 

These results are neither universal nor source- independent, but are useful generalizations of previous results to the 
case A\ ^ 4. 

The numerator of the first-order correction (|42[1 is identical to the left-hand side of Ij37|l , yielding 

9i = 0. (46) 

As with the total magnification, the first-order correction to the centroid vanishes universally in the PPN framework, 
which also means that centroid corrections beyond zeroth-order will be more challenging to observe directly. 

Remark. For the special case of the Schwarzschild metric in general relativity, Ebina et al. [25| and Lewis & Wang 
|20j found that the first-order corrections to the total magnification and centroid vanish. We have now generalized 
that result to all PPN models. 

4- Differential time delay 

In many cases the only observable time delay is the differential delay between the positive- and negative-parity 
images, 

Af = f--f+, (47) 

which we write as a series of the form 

Af = Af + Afie + 0{ef. (48) 
Starting from i|24H - 12{j [) and using (|27|) to write 9q in terms of (3, we obtain 

An = ^ (50) 

These relations are not universal in the sense that each depends on PPN parameters as well as the source position, 
but they are still useful for calculations. 
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IV. IMPLICATIONS FOR OBSERVATIONS 



In this section we employ the new lensing relations to assess prospects for using lensing observations to test theories 
of gravity. We consider three likely astrophysical scenarios: lensing by the Galactic black hole; conventional Galactic 
microlensing; and lensing in a binary pulsar system. 



A. Review of observables 



To facilitate this discussion, let us review the possible lensing observables, focusing on the weak-deflection limits 
plus the first-order corrections which should be measurable now or in the near future. To connect with realistic 
observations, we revert from our convenient mathematical variables (9, (3, /i, f) to true observable quantities B 7 
F, r). 

The traditional lensing observables are the positions, fluxes, and time delays of the images. The fluxes are related 
to the magnifications via the source flux: Fi = \^i\F src . (As an observable quantity, flux is positive-definite.) In 
principle, one could simply take the measured values, adopt the formulas derived in this paper as a model, and fit for 
the unknown parameters. However, with the relations found in Section IlIII as a guide, we believe it is instructive to 
make certain combinations of observables as follows: 



AF = F+ - F' 
tf+ F+ - 0- F~ 



F u 



At 



A^ 2 E + B 2 + £ + o (e) 2 . 



F tot = F+ + F- = F S1 



A^l + 2B 2 



— F — F 



2\B\^/A^ E +W 



+ 0(e)\ 



An 0? 



2(An?| + B 2 ) 3 / 2 



e + 0{e)\ 



= \B\ 



ZAx<d 2 E + AB 2 



2A X {) 2 + AB 2 



0{e) 2 



cd LS 



B\ X /A^ 2 E + B 2 



^\B\e + 0(ef 



(51) 
(52) 

(53) 

(54) 

(55) 



Ax 4% ln [ ^A 1 d 2 E + B 2 + \B\ 



^A 1 $ E + B 2 -\B\ 



(56) 



These equations summarize our results for PPN lensing, and play the key role in understanding how lensing can be 
used to test PPN theories of gravity. 



B. The Galactic black hole 

The center of our Galaxy is believed to host a supermassive black hole with a mass of M. = (3.6±0.2) x 10 6 Mq p^ : 
the distance to the lens is = 7.9 ±0.4 kpc |2^. Adopting the nominal values and neglecting the small uncertainties, 
we find the black hole's gravitational radius to be m. = 5.3 x 10 9 m = 1.7 x 10~ 7 pc, which corresponds to an angle 
of d, — 4.5 x 1CP 6 arc s. The corresponding lensing time scale is te = 71 s. 

We consider a source that is orbiting the black hole at a distance d^s <C d^ (so ds ~ d^). In the following 
quantitative estimates we let d* LS = rfis/(l pc) to simplify the notation. If the orbit is close to edge-on, part of it 
will lie close enough to the black hole (in projection) that the source can be significantly lensed. The angular Einstein 
radius is then $e — 0.022 (d Ls y^ 2 arc s, and our dimensionless perturbation parameter is e — 2.1 x 10~ 4 x (d^g)^ 1 ^ 2 . 

These numbers indicate that the two lensed images could be resolved with existing technology. From l|51|) . the 
angular separation of the two images is at least -J~A{ d E — 2 = 0.044 (d^g) 1 / 2 arc s (using Ai = 3.99966 ± 0.00090 
[23). At optical wavelengths, the Hubble Space Telescope has a resolution of about 0.05 arc s, while the CHARA 
interferometer |2^| can obtain a resolution of better than 10 -3 arc s. At radio wavelengths, interferometry can also 
achieve a resolution of ~ 10~ 3 arc s. The position uncertainties are much smaller than the resolution element; for 
example, radio interferometry observations of known lenses have yielded position uncertainties at the level of 10 -6 arc 
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s (e.g., [30| ) - With the conservative assumption that the images must be separated by 10 -2 arc s to be well resolved, 
we can still consider sources as close to the black hole as d^s = 0.05 parsec, and we can expect to measure the image 
positions with micro-arcsecond precision. 

A single observation could therefore be expected to yield the position and flux of each of the two images. Using 
equations (|5f fl - l|54[) . the four numbers (i9 + , $~ , F + , F~) would allow us to solve for four unknowns, which we could 
take to be cIls, B, F src , and A2. (Here we imagine taking the values of M., dr,, and A\ as given above.) Thus, in 
principle a single observation of an appropriate source could test gravity theories by measuring the invariant PPN 
bending-angle parameter A^. The position shifts that depend on A2 are of order $e £ = = 4.5 x 10~ 6 arc s, so 
existing technology has sufficient precision to measure them. Note that A2 is connected to the occurrence of naked 
singularities for certain gravity theories (see Sec. III.D of Paper I), so lensing could provide an observational test 
of the Cosmic Censorship conjecture [23. This analysis indicates that the main challenge for using lensing by the 
Galactic black hole to test gravity theories is just to find a source that is lensed. 

We could do even better with repeated observations, watching as the source moves and causes the images to change. 
We may estimate the time scale for such variations as the time it takes the source to move one linear Einstein radius, 

ds$E- For a circular orbit, this is Te — 6.5d* LS yr (independent of the black hole mass, since the Einstein radius 

1/2 

and the Keplerian orbital velocity both scale as Mj ). Keplerian orbital motion can be described with just five 
parameters: semimajor axis, period, eccentricity, inclination, and longitude of periastron. Repeated observations 
would thus allow us to determine a good model for 6 as a function of time. The universal relations then tell us 
that the quantity d + — $~ would be the most interesting combination of observables. From \5ty) . if d + — $~ has 
any dependence on B that is not strictly linear, that would represent a clear detection of higher-order effects from the 
gravity theory. The practical value of the universal relations, then, is to identify combinations of observables that 
would give the most direct evidence that the measurements are probing beyond the weak-deflection limit. 

The other reason to make repeated observations is of course to obtain more observables than unknowns, so the 
problem becomes overconstrained. In this case the data will not only determine the parameter values, but determine 
whether the PPN framework itself is an acceptable model. 



C. Galactic microlensing 

In conventional microlensing, a foreground star or compact object lenses a star in the Galactic bulge. In what follows 
we quote M. in units of the mass of the Sun, and d$ in units of 8 kpc: M* — M, / M© and d* s = ds/(8 kpc) . We consider 
a typical situation with the lens lying about halfway between the observer and source: d^ ~ d^s ~ d$/2. The angular 
gravitational radius is then i9. ~ 2.5 x 1CT 12 x (M*/d* s ) arc s, the angular Einstein radius is -d E ~ 10~ 3 (M*/^) 1 / 2 
arc s, and the perturbation parameter is e ~ 2.4 x 10~ 9 x (M*/^) 1 / 2 . 

Present microlensing programs only measure the total flux as a function of time (or equivalently source position). 
Equation (|53|l shows that there is no first-order correction to the total flux, so it is not feasible to test theories of 
gravity with microlensing at present. Future programs may be able to resolve the images (which will be separated 
by a few milli-arcseconds). But in order to test gravity theories they would need to measure image positions with a 
precision at the level of 10 -12 arc s, or fluxes with a fractional uncertainty of order e ~ 10~ 9 . In other words, it is 
not reasonable to expect to test theories of gravity with conventional microlensing in the foreseeable future. 



D. Pulsars in binary systems 

Hopes for using stellar-mass lenses to test theories of gravity are not lost. The amplitudes of the correction terms 
are governed by 



( GM. 



\ Ac 2 d L d L s 



1/2 

(57) 



so we may be able to use stellar-mass lenses if we can find systems where dLS is sufficiently small. The ideal system 
would be a pulsar in a binary system with a compact object (another pulsar, or a black hole), in an orbit seen nearly 
edg e-on. An example of such a system was recently discovered: the binary pulsar J0737— 3039 [3JJ. Rafikov & Lai 
|32| have made detailed calculations of various effects on the pulsar timing measurements, including not only multiple 
imaging but also relativistic aberration and latitudinal delays associated with the spin of the source. We use this 
system more generally to be representative of binary systems consisting of a pulsar and a compact object, and to 
illustrate the amplitude of lensing effects associated with different theories of gravity. 
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In J0737— 3039, we take the fast millisecond pulsar to be the light source, and the slow pulsar with M, = 1.25 Mq 
to be the lens. The binary orbit has a semimajor axis a — 8.78 x 10 5 km. The orbital eccentricity is fairly small 
(e = 0.088), so for illustration purposes take d^s = a. The lens gravitational radius is then m. = 1.8 km, so the lensing 
time scale is t e = 2.5 x 10~ 5 s. Using d L sa ds, the physical Einstein radius is (Il'&e ~ (4GAf.d is /c 2 ) 1 / 2 = 2.5 x 10 3 
km. (The angular Einstein radius cannot be determined because the distance to the lens is not known.) The 
perturbation parameter is e = 7.2 x 10~ 4 . 

The two images could not be resolved spatially. They could be resolved temporally, though: when the source is 
behind the lens, each radio pulse would actually consist of two pulses (one from each image) separated by the lens 
time delay (see 32]). The amplitudes of the two pulses could be measured; and the intrinsic pulse amplitude could be 
measured when the source is not behind the lens. Thus, in this system the observables would be F + , F~, and At as 
a function of source position B, as well as F src . Should one wish to go further, analysis of the source's orbital motion 
would yield a prediction for the pulse arrival time in the absence of lensing and make it possible to measure the time 
delays t ± for the two images separately. 

In this scenario, the universal relations indicate that for testing gravity theories the most valuable measurement 
would be the flux difference AF = F + — F~ as a function of source position. From l|54(l . this quantity is constant 
in the weak-deflection limit. Thus, any variation in AF with source position would reveal that the measurements are 
probing beyond the weak- deflection limit. Again we see the universal relations helping us identify combinations of 
observables that are best suited for testing gravity theories. 

V. CONCLUSIONS 

Using gravitational lensing by compact objects, we have presented new propects for testing theories of gravity 
within the PPN framework. In this paper we generalized the PPN lensing formalism from Paper I to include fully 
general third-order PPN models. We determined the weak-deflection limits plus first- and second-order corrections in 
e = t?,/i?£ for observable properties of lensed images (positions, magnifications, and time delays). 

During the PPN analysis, we discovered some surprising new fundamental relations between lensing observables. 
Some of the relations are universal for the entire family of PPN gravity models. A deep conceptual implication is 
that any observed violation of the universal lensing relations (given that assumptions [Al]-[A3] apply) would indicate 
that a fundamentally different theory of gravity is at work — one outside the PPN framework. The new relations 
have enabled us to identify combinations of lensing observables that are key to probing the spacetime metric by 
constraining the invariant PPN parameters of the light bending angle. The parameter Ai is related to the existence of 
naked singularities in certain gravity models (see Paper I), so constraining Ai also provides a possible observational 
test of the Cosmic Censorship conjecture. The new lensing relations will, in other words, play important roles in 
planning observing missions to test theories of gravity. 

In a practical application, we identified lensing observables that are accessible to current or near-future instrumen- 
tation, considering three likely lensing scenarios: the Galactic black hole, Galactic microlensing, and the binary pulsar 
J0737— 3039. A noted application of the new lensing relations is the ability to find combinations of observables that 
will yield a direct method for knowing when observations are probing beyond the standard weak-deflection regime. 
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